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Abstract 

For aerodynamic and hydrodynamic vehicles, it is highly desirable to reduce drag and noise levels. A reduction 
in drag leads to fuel savings. In particular for submersible vehicles, a decrease in noise levels inhibits detection. 
A suggested means to obtain these reduction goals is by delaying the transition from laminar to turbulent flow 
in external boundary layers. For hydrodynamic applications, a passive device which shows promise for transition 
delays is the compliant coating. In previous studies with a simple mechanical model representing the compliant wall, 
coatings were found that provided transition delays as predicted from the semi-empirical e n method. Those studies 
were concerned with the linear stage of transition where the instability of concern is referred to as the primary 
instability. For the flat-plate boundary layer, the Tollmien-Schlichting (TS) wave is the primary instability. In one 
of those studies, it was shown that three-dimensional (3-D) primary instabilities, or oblique waves, could dominate 
transition over the coatings considered. From the primary instability, the stretching and tilting of vorticity in the 
shear flow leads to a secondary instability mechanism. This has been theoretical described by Herbert based on 
Floquet theory. In the present study, Herbert’s theory is used to predict the development of secondary instabilities 
over isotropic and non-isotropic compliant walls. Since oblique waves may be dominant over compliant walls, a 
secondary theory extention is made to allow for these 3-D primary instabilities. The effect of variations in primary 
amplitude, spanwise wavenumber, and Reynolds number on the secondary instabilities are examined. As in the 
rigid wall case, over compliant walls the subharmonic mode of secondary instability dominates for low-amplitude 
primary disturbances. Both isotropic and non-isotropic compliant walls lead to reduced secondary growth rates 
compared to the rigid wall results. For high frequencies, the non-isotropic wall suppresses the amplification of the 
secondary instabilities, while instabilities over the isotropic wall may grow with an explosive rate similar to the rigid 
wall results. For the more important lower frequencies, both isotropic and non-isotropic compliant walls suppress 
the amplification of secondary instabilities compared to the rigid wall results. The twofold major discovery and 
demonstration of the present investigation are: (1) the use of passive devices, such as compliant walls, can lead to 
significant reductions in the secondary instability growth rates and amplification; (2) suppressing the primary growth 
rates and subsequent amplification enable delays in the growth of the explosive secondary instability mechanism. 


1. Introduction 

Research involving flow over flexible walls was 

started in the late-1950’s by Kramer 1,2 . Experimentally, 

Kramer found significant drag reductions using rubber 

coatings over rigid walls. Investigators in the 1960’s fo- 

cused on the task of experimentally duplicating and the- 

oretically explaining Kramer’s results. The majority of 

these studies failed to produce any comparable results; 

yet, the theoretical results laid the foundation for all fu- 

ture studies involving flexible walls. Interest turned to- 

ward the use of compliant walls for turbulent drag reduc- 
tion. In the 1970’s NASA 1 * 3 and in the 1980’s the Office of 

Naval Research 4 sponsored investigations involving the 
use of compliant walls for the turbulent problem. Al- 
thought most of the results from this era were either in- 

conclusive or unsatisfactory, the contributions, together 
with earlier results, have acted as stepping stones to the 
understanding of the physically 


complex fluid/wall interaction phenomena. A compre- 
hensive review of the pioneering studies was given by 
Bushnell, Hefner and Ash 5 , in particular for the tur- 
bulent flow problem. More recent reviews were given 
by Riley, Gad-el-Hak and Metcalfe 6 * * , Gad-el-Hak 7,8 , and 
Carpenter 9 * * . 

Motivation for the present investigation is partially 
derived from the following favorable theoretical and ex- 
perimental results. In the early 1980’s, Carpenter and 
Garrad 10,11 showed theoretically that Kramer-type sur- 
faces could lead to potential delays in transition. Fur- 
ther, they indicated deficiencies in previous investiga- 
tions which may have prevented their achieving results 
comparable to Kramer’s. Only recently, experiments 
performed by Willis 12 and Caster 13 showed favorable re- 
sults using compliant walls. As outlined in the above 
mentioned reviews, a number of investigations in the 
past ten years have been conducted involving flexible 
walls. A main emphasis of these studies was to under- 
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stand the physical mechanisms involved in the fluid/wall 
interaction of transitional and turbulent flows. Most of 
these studies focused on the 2-D instability problem, ex- 
cept Yeo 14 who showed that a lower critical Reynolds 
number existed for the isotropic compliant wall for 3-D 
instability waves. Carpenter and Morris 15 and Joslin, 
Morris and Carpenter 16 have shown that 3-D Tollmien- 
Schlichting waves can have greater growth rates over 
compliant walls than 2-D waves. However, they showed 
that, even though 3-D waves may be dominant, transi- 
tion delays are still obtainable through the use of com- 
pliant walls. They considered a compliant wall model 
used by Grosskreutz 17 for his turbulent boundary-layer 
experiments. 

In this paper, the growth rates and amplification 
of secondary instabilities over compliant walls are pre- 
dicted and compared to the rigid wall results. A the- 
ory for secondary instabilities is used which is based on 
Floquet Theory and was developed by Herbert 18 . This 
theory agrees remarkably well with experiments, in par- 
ticular those of Klebanoff, Tidstrom and Sargent 19 for 
peak-valley splitting (fundamental) and Kachanov and 
Levchenko 20 for peak- valley alignment (subharmonic). 
The secondary instability theory is extended to allow for 
3-D primary instabilities, which are dominant over the 
compliant walls considered. Appropriate boundary con- 
ditions for the Grosskreutz 17 wall model are derived for 
the secondary analysis. 

In the next section, the primary instability problem 
is outlined. In section three, the secondary instability 
problem is discussed. The numerical methods, results, 
and a summary follow in the remaining sections. 

2. Primary Wave Model 

Results for 2-D and 3-D primary instabilities over 
compliant and rigid walls have been presented by 
Joslin 21 , Morris and Carpenter 15,16 . Since the secondary 
instability theory is based on and includes a primary in- 
stability analysis, brief derivations of the dynamic equa- 
tions and boundary conditions for the 3-D primary insta- 
bility problem are included below. The disturbances are 
represented as travelling waves which may grow or decay 
as they propagate. Nonlinear coupling is ignored so that 
individual components of the frequency spectrum may 
be studied. Additionally, the quasi-parallel assumption 
is made. 

Consider an incompressible laminar, boundary-layer 
flow over a smooth flat wall. The Navier-Stokes equa- 
tions govern the flow. The Blasius profile is used to 
represent the mean flow. 

A small-amplitude disturbance is introduced into 
the laminar flow. A normal mode representation is given 
as 


M, fi'iK®, V, A I) ={vi, fii}(2/) exp[i(*a cos <f> 

+za sin <j> — uit)] + c.c. (1) 

where iq and fR are the complex eigenfunctions of nor- 
mal velocity and vorticity, respectively. To obtain a real 
solution, complex conjugate solutions denoted by c.c. are 
required, a is the wavenumber, u> is the frequency, and <f> 
is the wave angle. In general, a and u> are complex lead- 
ing to an ambiguity in the system. For temporal analy- 
ses, a is a real specified wavenumber and u> is the com- 
plex eigenvalue. For spatial analyses, u> is a real specified 
frequency and a is the complex eigenvalue. For the com- 
pliant wall problem, Joslin, Morris and Carpenter 16 have 
shown that the use of eqn. (1) leads to an overestimation 
of the growth of the wave as it propagates. The wave ac- 
tually propagates in a nearly streamwise direction which 
is in the direction of the group velocity, not normal to the 
wave fronts. In the present paper, the secondary insta- 
bilities are investigated using this simple representation 
of the primary instabilities. Since the present approach 
is conservative, it should exemplify the benefits of using 
compliant walls as a means to obtain transition delays. 
Also, a major emphasis and motivation of the present 
study is to determine the behavior, or response, of the 
phenomena-namely secondary instabilities-to compliant 
walls. 

If the normal mode relation (1) is substituted into 
the linearized form of the Navier-Stokes equations, the 
following nondimensional system results: 

v'l" + ai(y)v" + a 2 (y)vi = 0 (2) 

where 

a\(y) = — iRs(U 0 (y)a cos <f> — w) — 2a 2 
a 2 (y) = iRsa 2 (U 0 (y)a cos <f> - w) 

+iRga cos <f)U 0 "(y ) + a 4 

and 

+ a 3 (y)ft i + a 4 (y)v 1 = 0 (3) 

where 

a 3 (y) = - a 2 - iR s (U 0 (y)a cos <f> - w) 
a A (y) = - iR s a sin <f>U' 0 (y). 

The equations are nondimensionalized using the 
freestream velocity Coo, kinematic viscosity v, and 
an appropriate length scale. Convenient lengths for 
the boundary-layer scale with the ^-Reynolds number, 
R x = Ucyox/v. These include a thickness, 6, where 
the Reynolds number is defined R$ = R^J 2 and a 
boundary-layer displacement thickness, 6* , where R/,» = 
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17207 B}J~ . Eqns. (2) and (3) are referred to as the Orr- 
Sommerfeld and Squire equations, respectively. The sys- 
tem requires six boundary conditions. Requiring that 
the disturbance fluctuations vanish at infinity supplies 
three: 

n (y) , vi(y) , 0 as y — o o (4) 


For small displacements of an element out of equi- 
librium, the mechanical surface can be thought to move 
in a direction perpendicular to the rigid swivel-arm. The 
horizontal and vertical displacements (£1,771) are linked 
to the angular displacement (66) by 

£1 = f<5f? sin and 771 = £69 cos 9 (5) 


The remaining boundary conditions are determined from 
the compliant wall model. 

The compliant wall model used for the present pa- 
per was introduced by Grosskreut.z 1 ' in his experimental 
drag reduction studies with turbulent boundary layers. 
He suggested that the link between streamwise and nor- 
mal surface displacements would cause a negative pro- 
duction of turbulence near the wall. Although his results 
for the turbulent flow were disappointing, the surface 
does react to the fluid fluctuations in transitional flow in 
such a way as to reduce production of instability growth. 
Carpenter and Morris 22 have shown by an energy analy- 
sis how the many competing energy-transfer mechanisms 
are influenced by the compliant wall presence. Of note 
is the reduced energy production by the Reynolds stress 
which may cause the reduced growth rates. Further, 
Joslin, Morris and Carpenter 113 predicted that transition 
delays of 4-10 times the rigid wall transition Reynolds 
number were achievable with this coating. So the model 
has been extended to allow for a secondary instability 
analysis. 

The mechanical model consists of a thin, elastic 
plate supported by hinged and sprung rigid members in- 
clined to the horizontal and facing upstream at an angle, 
9, when in equilibrium. A sketch of the mechanical wall 
model is shown in Fig. 1. The boundary conditions are 
obtained by enforcing a balance of forces in the stream- 
wise and spanwise directions and the continuity of fluid 
and wall motion. These are given below in linearized 
form. 
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Fig. 1 Mechanical model representing the 
Grosskreutz compliant coating. 


where £ is the length of the rigid-arm member. Equations 
of motion for the element in the streamwise and spanwise 
directions may be obtained by a balance of the forces of 
the fluid fluctuations acting on the surface and the forces 
due to the wall motion. These equations are 


; d2 Vi , , n <9 4 ?7 i 
p m b df , +(5, Qxi 


dx 2 dz 2 


-) cos 2 9 + I\ E Vi - E x b 


5 2 £i . 


= (~P+ Tyy)< 


1 + T vx sin 9 cos ( 


,d 2 Ci , r . _ p , d 2 Ci 

' nr + ' nr 


where 




{E. r ,E : }b 3 

12(1 ^x^z) 


and B r ~ = 


£1 is the spanwise surface displacement. p m and b are the 
plate density and thickness; (S, F , 5. r; , B~ ) are the flexu- 
ral rigidities of the plate in the streamwise, transverse, 
and spanwise directions; (E x , E ~ ) are the moduli of elas- 
ticity of the plate; I\ E , Es are the effective streamwise 
and spanwise spring stiffness; p is the pressure fluctua- 
tion which is obtained from the fluid momentum equa- 
tions; and T yx , Tyy and T y ~ are the streamwise, normal, 
and spanwise viscous shear stress fluctuations in the fluid 
acting on the wall. 

The terms on the left hand side of eqn. (6) refer 
to mechanical forces and the terms on the right refer to 
fluid motion forces due to viscous stress and pressure 
fluctuations. For the case where the ribs are aligned 
at 6 = 0°, the wall becomes isotropic and reduces to the 
theoretical model studied by Carpenter and Garrad 10,11 . 
Otherwise the wall is referred to as non-isotropic and the 
rib angle is determined by 6. 

The continuity of fluid / wall motion is given in the 
streamwise, normal, and spanwise directions, respec- 
tively as 

^ = «i +ViK (8) 
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where («i, iq, uq) are the disturbance velocity compo- 
nents in the streamwise, normal and spanwise directions. 
For the Grosskreutz coating, Ks oo is assumed, which 
from eqn. (7) would result in zero effective spanwise 
surface displacement. This implies from eqn. (10) that 
uq(0) = 0. Strictly speaking, if the assumption Ks oo 
is relaxed, the resulting instabilities have larger growth 
rates. This suggests that spanwise stiffeners are stabiliz- 
ing to a disturbed flow. So with the assumption enforced, 
a better coating for potential transition delays results. 
The surface displacement takes the same normal mode 
form as the primary wave given by eqn. (1). The normal 
modes are substituted into eqns. (6-10). The equations 
can be reduced to three equations in terms of the normal 
velocity and vorticity 16,21 . 

3. Secondary Instability Theory 

In this section, a secondary instability theory devel- 
oped by Herbert 18 is extended to allow for 3-D primary 
instabilities. Additionally, boundary equations describ- 
ing the compliant walls are introduced for secondary in- 
stabilities. The flow is governed by the Navier-Stokes 
equations. Instantaneous velocity and pressure compo- 
nents are introduced and given as 

v(x, y, z, t) = v 2 (x, y, z, t) + Bv 3 (x, y, z, t) 

p(x, y, z, t) = p 2 (x, y, z, t) + Bp 3 (x, y, z, t) (11) 

where p 3 and v 3 = (u 3 , v 3 , w 3 ) are the secondary dis- 

turbance pressure and velocity in the fixed laboratory 
reference frame ( x,y,z ); and p 2 and v 2 = (u 2 , v 2 , w 2 ) 
are the basic pressure and velocity given by, 

v 2 (x, y, z,t) = {U 0 (y), 0, 0} + A{u x , iq, wq}(£, y, z,t) 
p 2 (x,y,z,t) = Api(x,y,z,t) (12) 

The basic flow is given by the Blasius profile and eigen- 
functions of the primary wave. Assume locally that the 
primary wave is periodic in t and periodic in (x, z) with 
wavelength A r = 27T j a r and ehne a disturbance phase 
velocity which is 

c r = (c x = ui r / a r cos <f>,0 , c z = ui r / a r sin <f>). 

Then in a frame moving with the primary wave, 

vi(x, y, z) = Ti(*, y, z) = Ui(x + \ x ,y, z + \ z ) (13) 

where (x, z) is the reference frame moving with the wave. 
With an appropriate normalization of primary eigen- 
functions («i,iq,uq) the amplitude, A, directly mea- 
sures the maximum streamwise rms fluctuation. This 
is given by 


The instantaneous velocities and pressure (11) are 
substituted into the Navier-Stokes equations which are 
linearized with respect to the secondary amplitude, B. 
The disturbance pressure is eliminated, resulting in the 
vorticity equations. 

[^ V2 - ^]^3-(T 2 • V)fi 3 - (v 3 ■ V)fi 2 

+(fi 2 • v)t 3 + (n 3 • V)v 2 = o (15) 

with the continuity equation 


V • v 3 = 0 


(16) 


As with the primary problem, the equations are formed 
in terms of the normal velocity (U 3 ) and vorticity (fi 3 ). 
The equations are found to take the form 21 , 


dU 0 dv 3 
dy dz 




+ ^{-(17! • V)fi 3 - (T 3 • V)fii 


dvi ,dv 3 dv 3 

+(tl+ ar>s 7 + fil ay 


dvi dv 3 
( Pl 8x > dz 


dv\ du 3 
dz dy 


dv\ dw 3 ' 
dx dy ■ 


= 0 


(17) 


and, 

r 1 „ 2 d d d 2 d 2 U 0 dv 3 

[_V=---(C„-c,) s + c.-]V^ + 




dy 2 dx 
d 2 ui d 2 ui 


d 2 vi d d 2 vi 
' Z - - n™ + I 


d 2 v\ d 2 v\ 


dz 2 

d 2 ui 


dy 2 


+ ( 


dxdy J dx ' v dz 2 dx 2 
d 2 w\ d 2 w\ d 2 u\ 


dy 2 


dx 2 


dxdz 
. dvi 


dy 2 
d 2 vi d 


dydz ’ dz + ^ dy 


)- 

dxdy dy 

dv\ 


d 2 


' dz 2 dx 2 dy 2 dx dxdy 
dv 


d 2 


_ 2 (— + — 1 - 

1 dz dx 'dxdz 


di ^ V3 + [ “4 (V2t;i) + 2( l^ 


2 d 2 ui d 
dxdy dx 

— ) + 2(^i 

dy 2 1 ^ dy 


2 ( 


d 2 wi 

dxdy 

2^1) 


d 2 vi d dvi d 2 

dxdz dz dx dx 2 


d 2 vi d 2 vi 
dx 2 dy 2 

-)2 Q2 

dz 2 


d 2 


dx 7 dxdy 


dz dGdJ M3+ t dz^^ + ^dydz 


} dwi d 2 
dx dydz 
d 2 ui d 2 v\ 


dvi d 2 
dz dx 2 


d 2 dui d 2 i \ 

dy 2 dz dxdy 3 J 


= 0 


)- 

dxdz 7 dx 

(18) 


max \ui{y)\ 2 

0<y <oo 


\ u i(ym)\ 


1/2 


(14) 
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The disturbance quantities i> 3 , Sl 3 and dv 3 /dy are re- 
quired to vanish far from the wall and at the wall for the 



rigid wall case. The compliant wall equations give the re- 
maining boundary conditions in the compliant case. Ad- 
ditionally, the primary amplitude, A, is a parameter in 
the equations and is assumed to be locally non-varying. 
As A —>■ 0 the Orr-Sommerfeld and Squire equations re- 
sult. For the case of interest where A zfz 0, the primary 
eigenfunctions (u i, tq, uq) appear in the equations as co- 
efficients. 

To solve the secondary problem, a normal mode so- 
lution having the following form is assumed, 

v 3 (x,y,z,t) = e at+i ^ zcos,i} - xsln ^V(x,y,z) (19) 


In the present study, temporally-growing tuned modes 
are examined. The temporal growth rate is <7 r , and oq 
can be interpreted as a shift in frequency. In this case, 
y r = ji = 0. If oq = 0, then the secondary disturbance 
is travelling synchronously with the basic flow. 

The boundary conditions for the secondary distur- 
bance are given as 

v n , v ' n , Cl n —> 0 as y oo (22) 

along with the compliant wall boundary conditions. In 
the rigid wall case, 


where (3 = 2ir/A z is a specified spanwise wavenumber 
and <7 = <j r T icq is a temporal eigenvalue or is spec- 
ified for spatial analyses. V(x,y,z) is a function that 
represents the class of secondary modes. Floquet theory 
suggests the form of solution for periodic systems. For 
the present problem, this may be written 

V(x,y,z) = e' l{xcos,t>+zsln ^V{x,y,z) (20) 


where 7 = y r + iji is the characteristic exponent and 
V(x, y, z) is periodic in the (x, z) plane and may be rep- 
resented by a Fourier series. Thus the representation of 
the secondary instability for a 3-D basic flow is, 


< 7 t-\-i(3(z cos (p — x sin cp)-\-j(x cos 4>-\-z sin <p~) 

c /3 — t 


Vniy ) 6 


i(n / 2)a r (x cos 4>-\-z sin <p~) 


(21) 


n = — 00 


This suggests a form of solution for the secondary distur- 
bance based on a coordinate system oriented at an angle 
<f> with respect to the mean flow and moving with the 
primary wave. If the coordinate system is aligned with 
the primary wave, or <f> = 0°, then the solution for the 
secondary disturbance would follow Herbert, Bertolotti 
and Santos 23 who considered a 2-D primary wave. 

If solutions given by eqn. (21) are substituted into 
eqns. (17) and (18), an infinite system of ordinary differ- 
ential equations result. The dynamic equations are de- 
termined by collecting terms in the governing equations 
with like exponentials. The system consists of two dis- 
tinct classes of solution because the even and odd modes 
decouple. Even modes correspond to the fundamental 
mode of secondary instability, and the odd modes are 
the subharmonic mode. Only a few terms of the Fourier 
series are retained since, as shown by Herbert, Bertolotti 
and Santos 23 , this provides a sufficiently accurate ap- 
proximation for a 2-D disturbance. 

This form of solution indicates two complex quanti- 
ties, it and 7, which leads to an ambiguity similar to 
that found with the Orr-Sommerfeld/Squire problem. 
There are four unknowns, a r , <7; , q r , 7; . Two can be de- 
termined while two must be chosen in some other way. 


v n , v ’ n , Cl n = 0 at y = 0 (23) 


The analysis for the compliant boundary conditions for 
secondary instabilities follows the same route as was 
taken for the primary instabilities, except a number of 
additional terms arise due to the presence of the primary 
wave. 

The fluid/wall motion must be continuous in each 
direction. In addition, the equations of force (6, 7) must 
balance in the streamwise and spanwise directions in the 
reference frame moving with the primary wave. Con- 
sistent with the fluid equations, the amplitude of the 
primary wave is assumed to be locally non-varying. In 
deriving the final form of the wall equations, a signifi- 
cant difference between the primary and secondary form 
arises from the pressure contribution. The pressure for 
the secondary disturbance is determined from the mo- 
mentum equations which are complicated by primary 
coupling terms. 

The continuity of motion between the fluid and solid 
is given by 


= « 3 + mK + a{(^ • v)« 3 + (€3 • v)«i} (24) 


_ w 3 -| - A-j (£ : • V)v 3 + (£3 

^- = w 3 + A < j (£i • V)w 3 + (£3 • 


where 


-r _ . d d d 

£ • V — £7: b y-z b C77" 

ox ay oz 


(25) 

(26) 


Equations (24-26) involve six unknowns for the velocity 
fluctuation and surface displacement in a highly coupled 
system. As with the primary boundary conditions, it 
is possible to derive a set of equations which represent 
the surface motion in terms of the normal velocity and 
vorticity only. This is algebraically very tedious. A com- 
plete derivation is given by Joslin 21 . Note that if A = 0 
in the secondary wall equations, the primary wall equa- 
tions result. This occurred with the fluid equations as 
well. 
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4. Numerical Methods of Solution 

The algebraic complexity of the dynamic equations 
for the secondary disturbance and the compliant wall 
equations requires that care be taken in applying any 
numerical technique. Because no theoretical or exper- 
imental data are available for the compliant problem, 
both shooting 16 and spectral 21 approximations are used. 
Also, noting that Bertolotti 24 has shown for the rigid 
wall problem with a 2-D primary instability that, after 
the transformation from spatial to temporal, the solu- 
tions are in good agreement, a temporal analysis is pre- 
sented in this paper. 

For the spectral method, Chebyshev series are intro- 
duced to approximate each mode of the Fourier series. 
An algebraic transformation is used to change from the 
Chebyshev spectral domain [—1, 1] to the physical do- 
main. Due to the properties of the Chebyshev polyno- 
mial, the equations are recast in integral form. Cheby- 
shev polynomials are used to represent the basic flow in 
the series which are substituted into the integral equa- 
tions. For the basic flow, 35 polynomials provided suffi- 
cient resolution of the eigenfunctions. The series repre- 
senting the secondary instability requires 40 polynomi- 
als for sufficient convergence to the dominant eigenvalue. 
For the shooting method, beginning with the equations 
for the compliant wall, integrations of the disturbance 
equations across the boundary layer are performed us- 
ing a Runge-Kutta scheme. At the edge of the bound- 
ary layer, the numerical solution vectors are matched 
with the asymptotic solutions. A very accurate initial 
guess is found to be required for convergence using this 
method. To demonstrate the accuracy of the numerical 
techniques, a comparison for the rigid wall case is made 
with Herbert 25 for Rg = 826.36, Fr = 83, j3 = 0.18, 
and A = 0.02. Herbert obtained the dominant mode 
<7 = 0.01184. In good agreement, the present spec- 
tral and shooting methods lead to a = 0.011825 and 
a = 0.011839, respectively. 

5. Results 

For all of the results that follow, the freestream ve- 
locity is 20 m/s, the density is 1000 kg/m 3 , and the kine- 
matic viscosity is 1 x 10 -6 m 2 /s. The coatings considered 
consist of both isotropic and non-isotropic walls. Both 
walls were optimized at R/,» = 2240 for 2-D primary 
instabilities. The isotropic wall has properties 9 = 0°, 
b = 0.735mm, E x = 1.385MN/m 2 , K = 0.354GN/m 3 
and p m = 1000 kg/m 3 ; and the non-isotropic wall has 
properties 9 = 60°, b = 0.111mm, E x = 0.509MN/m 2 , 
K = 0.059GN/m 3 and p m = 1000 kg/m 3 . A Reynolds 
number of 2240 was chosen because, for a boundary layer 
over a rigid wall, the disturbance with the critical fre- 
quency (in the e n sense) reaches its maximum growth 
rate near this value of Reynolds number. Accordingly, 


this is a good choice of Reynolds number for optimizing 
the wall properties. In considering 3-D instabilities, the 
walls optimized for 2-D instabilities are used with the ad- 
dition of isotropic plates. The properties of an isotropic 
plate are direction independent; that is, E x = E z . Al- 
though complete details of the optimization process and 
philosophy are given by Carpenter and Morris 22 , a recap 
follows. 

With a flexible wall present, other modes of instabil- 
ity arise. With changes in the compliant wall properties, 
stable, or marginally stable, fluid and wall modes can be- 
come unstable and dominant. The present wall proper- 
ties were varied to achieve an optimal specified condition. 
This desired condition was to achieve a minimum growth 
rate for a dominant 2-D Tollmien-Schlichting instability 
while keeping other modes marginally stable. For the 
secondary analysis, these “optimal” compliant walls led 
to no additional unstable modes. However, this is not to 
say that additional growing modes may not appear for 
different wall properties. 

In this analysis, the primary wave amplitude (A) 
and the secondary instability spanwise wavenumber (/?) 
are parameters of the problem. Herbert 26 showed in his 
boundary-layer studies for the 2-D primary wave over a 
rigid wall that as the amplitude increases the growth rate 
of the secondary instability increases. Also, as the span- 
wise wavenumber is varied, the temporal growth rate 
reaches a maximum for a particular wavenumber. Ad- 
ditionally, Herbert showed that the subharmonic mode 
reaches greater growth rates than the fundamental mode 
for low amplitude disturbances. These findings were 
verified by the direct numerical simulations of Spalart 
and Yang 27 . Althought both the subharmonic and fun- 
damental modes over the compliant walls were exam- 
ined, emphasis is placed on the subharmonic mode, since 
as both theory and computations indicate, subharmonic 
disturbances are more unstable than fundamental distur- 
bances for small amplitudes. Limited fundamental dis- 
turbance results are included to verify that these modes 
do not become the dominant instability over compliant 
walls. 

Primary waves with frequencies that give maxi- 
mum disturbance growth rates are considered. For the 
isotropic wall, the maximum growth rate occurs at a fre- 
quency us = 0.065 (Fr =~ 29.0), where Fr = us/Rx 10 6 . 
Figure 2 shows the growth rates of the subharmonic and 
fundamental disturbances as a function of the spanwise 
wavenumber for the rigid wall and isotropic compliant 
wall. As the figure shows, growth rates over the com- 
pliant wall are reduced in comparison with the rigid 
wall results over the whole range of spanwise wavenum- 
bers. Additionally, the subharmonic disturbance has 
much larger growth rates than the fundamental distur- 
bance, as expected. Similar trends are found in the com- 
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parison of non-isotropic and rigid wall results. In consid- 
ering reductions in the growth rates of the subharmonic 
mode as a result of compliant walls, the isotropic wall 
suppressed the maximum growth rate by 20%). For the 
maximum growth rate over the non-isotropic wall, which 
occurs at a difference frequency than the isotropic wall 
case, the non-isotropic wall led to a reduction of 17%) 
compared to rigid wall results. So for a fixed Reynolds 
number and primary wave amplitude, both isotropic and 
non-isotropic compliant walls lead to reduced secondary 
instability growth rates compared to the rigid wall re- 
sults. 



Fig. 2 Growth rates of the secondary instabil- 
ities as a function of spanwise wavenumber for 
Rf = 2240, Fr ~ 29.0, and A = 0.01. subharmonic: 
— o— , rigid wall; • • o- • -, isotropic wall and funda- 
mental: — x — , rigid wall; • • • x • • •, isotropic wall. 

A more revealing measure of the effectiveness of us- 
ing compliant walls to suppress secondary instabilities is 
to compute the amplitude growth and decay with down- 
stream distance. The amplification of the primary and 
secondary instabilities are governed by 

In = [ —a*dx and In — = [ —dx (27) 

where A 0 ,B 0 are the initial amplitudes at x 0 , A and B 
are the amplitudes at a downstream distance x, and * 
denotes dimensional quantities. As shown by Herbert. 18 , 
the theoretical prediction of primary and secondary am- 
plification by eqns. (27) compares well with the exper- 
iments of Kachanov and Levchenko 20 . A similar com- 
parison is shown in Fig. 3, where the theoretical results 
were obtained with the present numerical techniques. 

Before computing similar amplification results over 
the compliant walls, inferences of the secondary insta- 
bility growth may be drawn from primary instability 
results. Figure 4 shows the maximum amplification of 
various frequency primary waves propagating over the 
rigid and compliant walls along with the amplification 


of waves at a frequency Fr = 53. At this frequency, 
primary amplitudes over the the isotropic wall are sim- 
ilar to those over the rigid wall. This suggests that the 
development of secondary instabilities might also be sim- 
ilar, since secondary disturbances are parametrically de- 
pendent on the basic flow. For the non-isotropic wall 
results shown in Figure 4, the primary amplitudes are 
suppressed significantly compared to the rigid wall and 
isotropic wall cases. Therefore, the non-isotropic wall 
would likely lead to a very different secondary instability 
development, most probably with reduced amplitudes. 
Again from Figure 4, one might expect greater differ- 
ences in the secondary instability development over both 
compliant walls as the Reynolds numbers increase and 
corresponding frequencies decrease. 
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Fig. 3 Amplitude growth as a function of 
Reynolds number for the subharmonic mode (B) 
of a 2-D primary wave (A) over a rigid wall at 

Fr = 124, A 0 = 0.0044, B 0 = 1.86x10-®, and b = 0.33. 
, theory and (x,o), Kachanov and Levchenko 

[15]- 



Fig. 4 2-D curves of maximum amplification for 

TSI waves over a , rigid wall; , isotropic 

wall; and , non-isotropic wall and — • — , waves 

of Fr ~ 53. 
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To demonstrate these postulations for secondary in- 
stabilities, eqns. (27) are used to compute the amplifi- 
cation of primary and secondary instabilities over the 
rigid, isotropic, and non-isotropic walls for A 0 = 0.004, 
B 0 = 1 x 1 0 — 5 , and b = 0.15, where b = j3/ R x 10 3 . The 
initial amplitudes are somewhat arbitrary: the present 
values were selected to be closgi- to the experiments of 
Kachanov and Levchenko 20 for rigid walls. The spanwise 
wavenumber (b) was chosen near the maximum growth 
rate of the secondary instability at the branches of the 
neutral curve for the rigid wall case. Both primary and 
secondary amplifications are shown in Fig. 5. Clearly, 
the results of the secondary instability growth over the 
isotropic and rigid walls are similar, as postulated, while 
the non-isotropic wall significantly suppresses the sec- 
ondary instability growth. 
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Fig. 5 Amplitude growth as a function of 
Reynolds number for the subharmonic mode (B) 
of a 2-D primary wave (A) at Fr ~ 53, A 0 = 0.004, 

B 0 = 1 x 10 -5 , and b = 0.15 over a , rigid wall; 

, isotropic wall; and , non-isotropic wall. 

The growth of the secondary instability is dependent 
on the parameters of the basic flow, most probably the 
primary instability amplitude (A). For example, Fig. 6 
shows the amplification of the primary and secondary 
instabilities over the rigid wall with properties as before 
and over the non-isotropic wall with both A 0 = 0.004 
and A 0 = 0.008. Even by doubling the initial amplitude 
of the primary disturbance, the growth of the secondary 
instability over the non-isotropic wall continues to be 
suppressed and has not exceeded the primary amplitude 
upon crossing the neutral curve. Yet,, the doubled initial 
amplitude (A 0 = 0.008) results in a significant increase 
in the secondary instability growth compared with the 
lower amplitude (A 0 = 0.004) results. Hence, the sup- 
pression of the primary instability amplitude is of utmost 
importance to suppress the onset of the secondary insta- 
bility growth. 
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Fig. 6 Amplitude growth as a function of 
Reynolds number for the subharmonic mode (B) 
of a 2-D primary wave (A) at Fr ~ 53, B 0 = 

1 x 10“ 5 , and b = 0.15 over a , rigid wall 

with A 0 = 0.004; — • — , non-isotropic wall with 

A 0 = 0.004; and , non-isotropic wall with 

A 0 = 0.008. 

.006 


£.004 

! 

1 

002 


000 

0.0 0 . 2 O A 0.6 

Tiivnumr^r. p 

Fig. 7 Growth rates of the subharmonic distur- 
bance for 2-D and 3-D primary waves over the 
isotropic wall as a function of spanwise wavenum- 
ber for R /,* = 2240, Fr ~ 22.3, and A = 0.01. , 

<j> = 0°|> ••,</> = 10°; and , <f> = 20°. 

Proceeding with investigating the effect of compli- 
ant. walls on secondary instabilities, 3-D primary waves 
are introduced and are determined by the specified wave 
angle In Fig. 7, subharmonic disturbance growth 
rates over the isotropic compliant wall are shown with 
variation in spanwise wavenumber and primary wave an- 
gle A frequency u> = 0.05 (Fr ~ 22.3) is selected 
since the 3-D primary wave growth rate is maximized 
at this frequency for the Reynolds number Rf = 2240. 
Subharnronic growth rates arising from the 2-D primary 
wave (<ft = 0°) are clearly larger than those from the 
3-D waves. As the primary wave angle (<ft) increases, 
the subharnronic growth rates continually become more 
damped. Additionally, for oblique waves (<ft yl 0°) the 
secondary disturbances no longer travel synchronously 
with the primary wave. This is shown in Fig. 8 by the fre- 
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quency shifts that result over the isotropic wall. Similar 
results occur for the non-isotropic wall. It is likely that 
this shift leads to a reduced efficiency of energy transfer 
from the basic flow to the secondary disturbance. This 
frequency shift is as much a result of the 3-D nature of 
the basic flow as it is of the compliant wall influence. Re- 
sults similar to those of Fig. 8 were found by Balachan- 
dar, Streett and Malik 28 for the rotating-disk problem 
which also has a 3-D basic flow. 
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Fig. 8 Frequency shift of the subharmonic distur- 
bance for 3-D primary waves over the isotropic 
wall as a function of spanwise wavenumber for 

R s , = 2240, 22.3, and ,4 = 0.01. - - - , 0 = 10° 

and , <f> = 20°. 


The secondary growth rate comparison in Fig. 7 is 
misleading as a result of holding the primary amplitude 
(A) fixed. With a variation in wave angle (<^>), corre- 
sponding primary wave amplitudes result as shown by 
Joslin, Morris and Carpenter 18 . Then a more realistic 
comparison of secondary growth rates arising from 2-D 
and 3-D primary waves should involve amplitudes suited 
to a given wave angle (<f>). One means to derive these 
amplitudes is through the use of eqns. (27). In Fig. 9, 
a comparison of the subharmonic growth rates over the 
isotropic compliant wall is made for 2-D and 3-D primary 
waves allowing for amplitude differences. Similar results 
are found for the non-isotropic wall as well. As shown, 
the 3-D primary wave does lead to much larger growth 
rates than a 2-D wave for the compliant walls considered. 
This clearly demonstrates that in spite of the loss of syn- 
chronization with the basic flow, 3-D primary waves lead 
to dominant secondary instabilities over compliant walls. 
Again, the amplitude of the primary wave is of utmost 
importance for determining the secondary disturbance 
growth. Yet, compliant walls do lead to a reduction in 
secondary growth rates compared with those for the rigid 
wall. Also, recall that the amplitudes for 3-D primary 
waves used in these calculations were determined by a 
normal mode assumption which leads to a conservative 


estimation of the primary amplitude lb , and a more re- 
alistic lower amplitude would lead to secondary growth 
rates somewhere between the 2-D and 3-D results shown 
in Fig. 9. But, the goal here is to determine the fun- 
damental effect of compliant walls on secondary distur- 
bances. 

Final amplitude calculations using eqns. (27) are 
carried out for the frequency Fr ~ 53. In Fig. 10, sec- 
ondary amplitudes arising from the most amplified 3- 
D primary wave over a non-isotropic compliant wall are 
compared to the previous results from 2-D primary waves 
over rigid and non-isotropic walls. At this frequency, the 
amplitude difference between the 2-D and 3-D primary 
instabilites is small, yet the secondary disturbance re- 
sponds notably. This is an indication that small changes 
in the primary instability, however slight, have a mount- 
ing effect on the rapidly developing secondary instability. 
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Fig. 9 Growth rates of the subharmonic distur- 
bance for 2-D and 3-D primary waves over the 
isotropic wall as a function of spanwise wavenum- 
ber for Rf = 1760 and Fr ~ 30.2 for , <f> = 0° 

with A = 0.010; •••,</> = 45° with A = 0.031. 
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Fig. 10 Amplitude growth as a function of 
Reynolds number for the subharmonic mode (B) 
of a 2-D and 3-D primary wave (A) at Fr ~ 53, 

B 0 = 1 x 10 -5 , A 0 = 0.004, and b = 0.15 over a , 

rigid wall — • — , 2-D non-isotropic wall; and • • •, 
3-D non-isotropic wall. 
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6. Summary 

In earlier studies 15,16,21 , it was shown that 3-D pri- 
mary instabilities theoretically dominate transition over 
the compliant walls considered, yet transition delays 
were found compared to the rigid wall. The present pa- 
per has further extended the understanding of the effect 
of compliant walls on transition mechanisms in bound- 
ary layers. Namely, the effect compliant walls have on 
secondary instabilities has been investigated. It has been 
shown that the use of compliant walls can lead to reduced 
growth rates and amplification of secondary instabili- 
ties. From both the earlier 3-D primary results combined 
with the present secondary findings, it should be empha- 
sized that the physical nature and make-up of the mech- 
anisms in transition are not altered by the control device 
(i.e. compliant wall). Rather, only the response of that 
mechanism is changed. This fact is of particular impor- 
tance for designing Laminar Flow Control (LFC) studies. 
As an example, the behavior of a secondary instability 
growth with variation in primary amplitude is well doc- 
umented by Herbert. As the primary amplitudes are 
reduced, the excitement of the secondary instability is 
delayed. Thus, active or passive devices which suppress 
primary instability growth should lead to corresponding 
suppression and delay of succeeding instabilities. This 
has been demonstrated above with the compliant wall. 
The twofold major discovery and demonstration of the 
present investigation is: (1) the use of passive devices, 
such as compliant walls, lead to significant reductions in 
the secondary instability growth rates and amplification; 
(2) suppressing the primary growth rates and subsequent 
amplification enable delays in the growth of the explosive 
secondary instability mechanism. 
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